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Abstract: Giant magnon solutions play an important role in various aspects of the 
AdS/CFT correspondence. We apply the dressing method to construct an explicit 
formula for scattering states of an arbitrary number of magnons on M x 5*^. The 
solution can be written in Hirota form and in terms of determinants oi N x N matrices. 
Such a representation may prove useful for the construction of an effective particle 
Hamiltonian describing magnon dynamics. 
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1. Introduction 

Classical string solutions in AdS^^ x play an important role in understanding vari- 
ous aspects of the AdS/CFT correspondence (see [1] for review). Integrability [2] is a 
powerful computational tool which has enabled many quantitative checks of the cor- 
respondence. A lot of work has been done exploring both string theory and gauge 
theory sides of the correspondence, culminating in the proposal for an exact S-matrix 
for planar M — ^ Yang-Mills theory [3] . 

Magnons are building blocks of the spectrum in the spin chain description of 
AdS/CFT. The Hofman-Maldacena elementary magnon corresponds to a particular 
string configuration moving on an M x S'^ subspace of AdS^ x S'^ [4]. String theory 
on M X 5^ (or M x 5*^) is classically equivalent to sine-Gordon theory (or complex 
sine-Gordon theory) via Pohlmeyer reduction [5, 6] (see [7] for AdS case). Giant one- 
magnon solutions on R x S*^ and M x S*^ map to one-soliton solutions in sine-Gordon 
and complex sine-Gordon respectively [4, 8]. Using this map, the scattering phase 
of two magnons was computed in [4] and shown to match that of [9]. Moreover, a 
sine- Gordon-like action has been proposed for the full Green-Schwarz superstring on 



AdS^ X [10, 11]. 
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In sine-Gordon theory, the dynamics of A^-sohtons is captured by the Ruijsenaars- 
Schneider model [12, 13]. Specifically, the eigenvalues of a particular N x N matrix 
entering into the description of the A'"-soliton solution (or r-function) of sine-Gordon 
evolve according to the Ruijsenaars-Schneider Hamiltonian. Positions and momenta 
in the Hamiltonian are related to the positions and rapidities of the solitons, and the 
phase shift for soliton scattering can be calculated from the quantum mechanical model. 
It is natural to wonder what the analagous Hamiltonian in the case of complex sine- 
Gordon and giant magnons is. Explicit A^-soliton solutions (in r-function form) serve 
as a useful starting point in deriving the Ruijsenaars-Schneider model from the sine- 
Gordon theory, and it is likely that a similar technique may prove useful for complex 
sine-Gordon and giant magnons as well. 

Interest for an effective particle description of giant magnon scattering emerged 
through the work of Dorey, Hofman and Maldacena [14] , where they illuminated the 
nature of double poles appearing in the proposed S-matrix of planar A/" = 4 Yang- Mills 
[3]. They were able to interpret these double poles as occurring from the exchange of 
pairs of particles, and in particular to precisely match their position on the complex 
domain with the prediction of [3], under the assumption that the exchanged particles 
are BPS magnon boundstates [15]. By studying the quantum mechanical problem 
corresponding to an effective particle Hamiltonian describing the scattering of two 
magnons with very small relative velocity, one should obtain an S-matrix whose double 
poles compare to the aforementioned results in the appropriate limit. 

Superposing magnons is a difficult problem because of the nonlinear equations of 
motion they satisfy. Integrability allows the use of algebraic methods such as dressing 
to construct solutions of nonlinear equations of motion [16, 17]. Indeed, the dressing 
method was used to describe the scattering of two magnons and spikes on M x S*^ (and 
various subsectors) as well as spikes in AdS^ [18]-[22]. However, it is a tedious process 
to obtain even the three-magnon solution. In this paper we will present an explicit 
string solution on R x 5"^ describing scattering of an arbitrary number of magnons 
by solving the recursive formula following from the dressing the {N — l)-magnon. 

The paper is organized as follows. In section 2 we review the dressing method for 
M-X and derive a recursive formula for the A/^-magnon solution in terms of {N — 1)- 
magnons. In section 3 we solve this recursion and present the A^-magnon solution. 
The solution can be presented in various ways, we find useful Hirota and determinental 
forms. As a consistency check we verify that our solution separates asymptotically into 
a linear sum of A^ well-separated single magnon solutions and demonstrate that the 
only nontrivial effect of the A^-magnon interaction is the expected sum of two-magnon 
time delays. The appendix clarifies the rules to construct the A'"-magnon solution and 
some examples are presented. 
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2. Giant magnons on M x 5^ 



The classical action for bosonic strings on IR x 5"^ can be written as 



^jdtdx [d'^X^'daX^ + A{Xi ■ X, - 1)] , (2.1 



^--2 

where ji runs from to 4 and i from 1 to 4. The Xi are embedding coordinates on 
and the Lagrange multiplier A constrains them on S^. 

After we impose the gauge X^{t,x) — t, eliminate A in terms of the embedding 
coordinates and switch to light-cone worldsheet coordinates z — {x—t)/2, z — {x+t)/2, 
the equations of motion and Virasoro constraints become 

ddZi + ^{dZjdZj + dZjdZj)Zi = 0, ZiZi = 1, (2.2) 

and 

dZidZi = BZidZi = 1, (2.3) 
where we have used the parametrization 

Zi = + iX2, Z2^Xs + iX^. (2.4) 

Giant magnons on IR x 5*"^ arc defined as solutions to the above system of equations, 
obeying the boundary conditions 

Zi{t,x^±oo) = e^*±^f/2+ia^ 

Z2{t,x^±oo) ^0. (2.5) 

The physical meaning of the boundary conditions is that the endpoints of the string lie 
on the equator of the on the Zi plane moving at the speed of light, and the quantity 
p called total momentum represents the angular distance between them. Finally, a can 
be any real constant. 

2.1 Review of the dressing method 

The dressing method is a general procedure for constructing soliton solutions to in- 
tegrable differential equations first developed by Zakharov and Mikhailov [16, 17]. It 
was applied in the context of giant magnons by some of the authors [18, 19], providing 
classical solutions for a variety of backgrounds. In what follows, we will review the 
basic steps of the method as they apply to the particular case of R x 5"^. 
We start by defining the matrix- valued field 

g{z,z)^(f-^-i^AeSU{2) (2.6) 
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and recasting (2.2) into 

dA + dB^ 0, (2.7) 
where the currents A and B arc given by 

A^zdgg-\ B^idgg-\ (2.8) 

The Virasoro constraints (2.3) can be also written as 

TrA^ = TrB^ = 2. (2.9) 

The nonhncar second order equation for g in (2.7) is equivalent to a system of 
linear first order equations for auxiliary field ^{z, A) 

A'^ - B^ 
id^ = ^ idyif = (2.10) 

1-A' 1 + A ^ ' 

provided (2.10) holds for any value of the new complex variable A called the spectral 
parameter, with A and B independent of A. 

Given any known solution g, we can determine A, B and solve (2.10) to find \1'(A) 
subject to the condition 

*(A = 0)=^. (2.11) 

Any ambiguity on factors that don't depend on z, z is removed by also imposing the 
unitarity condition 

[*(A)]^*(A) = 7. (2.12) 

It is easy to show that the equations of motion for the auxiliary field (2.10) are covariant 
under the following transformation with a A-dependent parameter 

*(A) ^ *'(A) = x*(A), 
A ^ A' ^xAx''+i('i- + ^)dxX-\ (2-13) 
B ^ B' =x5x-'+^(l-A)axx-^ 

under the condition that A' , B' remain independent of A. Thus, performing the above 
transformation to the known solution {^{\),A,B) produces a new solution to (2.7) 
with g' = ^'(A = 0). 

The condition (2.12) imphes that x(A) must obey 

[X(A)]^X(A) = /, (2.14) 

whereas the demand that A', B' are independent of A can be translated as further 
constraints on the analytic properties of x(A)- For the R x 5"^ case it turns out [18] 
that the dressing factor x(A) is 

x(A) = / + yf^P, (2.15) 
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where Ai is an arbitrary complex number and the hermitian projection operator P is 
given by 

P=^, i;i = *(Ai)e, (2.16) 

where e is an arbitrary vector with constant complex entries called the polarization 
vector. The projector P does not depend on the length of the e vector. 
The determinant of x(A) is 

det x(A) = (2.17) 

and if we want our dressed solution x(0)^'(0) to sit in SU{2) we should rescale it by 
the compensating factor \/\if\i. 

Putting everything together, the new solution g' — ^'(A = 0) to the system (2.7) 
is given by 




i(/ + :^p),. (2.18) 



2.2 Application and recursion 



This procedure can be repeated with g' as the solution we begin with, in order to obtain 
another new solution. In fact, once we have solved the differential equation (2.10) for 
^'(A) the first time, we no longer need to repeat this step for ^'(A), as we have that 
information already. Thus, from this point the method proceeds iteratively in a purely 
algebraic manner. 

More specifically, we can show that the auxiliary field ^^(A) that is constructed 
after N iterations is related to the auxiliary field ^^~^(A) occuring after A^— 1 iterations 
through 



^ ^1 A^ (A - X^){ab - cd) ^l^g ' ' ^'-'^^ 



where 



< = {-cd{X - Xn) + ab{X - AAr)Xr'(A) - ac{XN - Xn)K''\>^), 
^^^2 = i-cdiX - Ajv) + ab{X - Ajv))*f2"'(A) - ac{Xr, - XN)^g-\X), 

< = {ab{X - Xn) - cd{X - XN)Mf\X) + KAiv - Aiv)*fr'(A), 

iPg = {ab{X - Xn) - cd{X - Ajv))*^2~'(A) + K^jv - Ajv)*f2"'(A), (2.20) 
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and 



d^K-\\^) + ^g-\X^). (2.21) 



The new solution of (2.7) follows from (2.19) when taking A = 0. Due to (2.6) we 
can then read off the relation between the coordinates of the two solutions as 

\Am\{(io — cd) 

^2 = mA 1^ [(«^^^ - cdXN){-tZt') + «c(A^ - Xn)Z^-'] . (2.22) 

|Ajv| [ab — cd) 

Starting with the simple 'vacuum' solution representing a point particle rotating 
around the equator in the Zi plane, 

Z2 = 0, (2.23) 

and using the polarization vector e = (1, 1) the dressing method yields [18] the single 
magnon solution on M x 5"^ first obtained in [8] as a generalization of the original 
Hofman-Maldacena giant magnon solution on M x 5*^. Applying the method once more 
using the same polarization vector as before then gives a solution which asymptotically 
reduces to a sum of two single magnon solutions, and whose conserved charges arc sums 
of the respective charges of two single magnon solutions. Hence it can be interpreted 
as a scattering state of two single magnons. 

Prom the above considerations, it is natural to expect that the A^-times dressed 
solution will correspond to a scattering state of N magnons. The quantities Aj are 
parameters of the A"-magnon solution which we can more conventionally express as 
Xi — rjc'^'/^, with Pi the momentum of each constituent magnon and n a quantity 
associated to its U{1) charge. 



3. The A/^-magnon solution 

Successive apphcation of the dressing method suggests a compact closed form for the 
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A'^-magnon solution, which can be written as follows 



^1 



Z-y = -i 



nf=i \H D ' 



(3.1) 



with 



D = exp 

/ii=0,l 



2N 



2N 



where 



Mi =0,1 

w=o,i 



i<j i=l 
2N 2N 

Y BijiniHj + {Hi - - 1)] + XI l^iC^iZi + Ci) 

i<j 1=1 
2N 2N 

X + {fii - - 1)] + '^[(Xi{2iZi) + {(li - l)Ci 



(3.2) 



i<j 



1=1 



Z 



+ 



z 



A, + l' 

e^'^ = - \j, 



(3.3) 



and iV is the number of magnons. 

In the above formula the indices i^j take the 2A^ values (1, 1, 2, 2, TV), i < j im- 
plies this particular ordering, and we identify A;^ = A^, Z^ = Z}^^ . The symbol X];Uj=o 1 
implies the summation over all possible combinations of //i = 0, 1, /xi = 0, 1, . . . , //jv = 
0, 1 under the conditions 



2N 



for iVi, D, 
iV+L foriVa. 



(3.4) 



This description makes contact with a variety of A'"-soliton expressions of other inte- 
grable systems (for example see [23]). 

We have numerically checked (3.2) for high number of magnons, whereas in fig. 1 
we plot 1 for the first 4 magnons. In appendix A we give some examples. 

Our M X 5"^ A^-magnon solution is reduced to the M x 5"^ one if we let the spectral 
parameters A; lie on a unit circle, — 1. 



^Alternatively we may define new quantities pk such that pii-\ = Xi and p2i = Xi, and similarly 
for Zi. These will take values 1, 2...2N as usual. 
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3.1 Hirota form of the solution 



It is possible to write Zi, of (3.1) in an equivalent form similar to Hirota's [24], 
where A^2, -D are given by 



-D = ^ d{%x, ^2, • • • , In) exp \li{Zi^ + H V ^i^)] , 

2nCn 

Ni = Y ni{ii,i2,...,iN)exp[2i{Zi^+Zi^^ ^^ijv)]) (3-5) 

2JvC'jV 

N2= ^ ^2(21,^2, • • • , iN+i) exp [2i{Zi^ + H h ^i^+J] , 

2JvC'jv+l 



- 8 - 



and 

(N) (AT) 

d{ii,i2,...,iN) = Yl n 

k<l<N N<m<n 

N (N) (N) 

ni{ii,i2,...,iN) ^YIK, n n ^^-^^ 

j=l k<l<N N<m<n 

2N (N+l) (N-l) 

n2{ii,i2,---,iN+i) ^ Yl \j n ^^^^i n 

j=Ar+l fe<KAr+l Ar+l<m<n 

where N is the number of magnons, indicates summation over all possible combi- 
nations of n elements taken from Yt^^ indicates the product of all possible combina- 
tions of the n elements, and \ij — Xi — Xj. Finally, we have arranged our 2N elements 
Zi as {Zi, Zi, . . . , Z]^} and our 2N A's as {Ai, Ai, . . . , X^}. We always assume that 
ii < . . . < iN- 

Finally, we should mention that we can get a more symmetric yet complicated- 
looking version of our A^-magnon expressions, by factoring out the terms 

nil Xi exp hi Xlii from A^i, 

< nil A/ exp (2^ Eii 2i) from N2, (3.7) 

exp X^ili from D. 

Written in this way, D has the nice feature of being real. More importantly, and 
as we will see in the following sections, this form of the A^-magnon solution is useful 
for analyzing its asymptotic behavior and demonstrates the symmetry that will allow 
us to write it in a determinant form. 

3.2 Determinant form for Zi 

It is known that for the (complex) sine-Gordon equation and several other integrable 
equations, the A'"-soliton expressions similar to (3.1)-(3.4) and (3.5)-(3.6) can also be 
rewritten in a form involving determinants of N x N matrices [25]. It is precisely 
expressions of this type that become particularly useful when extracting the effective 
particle description of the soliton problem [13]. Motivated by the same goal for the 
case of giant magnons, we haven been able to find a determinant formula for Zi. In 
particular, we may write 

— ^XiV^^ det{I + A-^FAF) 



det(7 + FF) ' ^^-^^ 
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where A, F are N x N matrices^ with elements 



Aw 


— SkiXi, 


Fki 


= e-^'^>'G, 


Gkl 


"TT Xkfh 



kl, 



(3.9) 



/c, / = 1, 2, . . . , A^, and / the identity matrix. Interestingly, the matrix G can further 
be expressed as G — H (ff) ^ where if is a matrix with elements H^i = (A^)'"^. The 
determinant of H is what is known in the literature as the Vandermonde determinant, 
given by the simple formula 

detH^HiXi-Xk). (3.10) 

k<l 

This decomposition in terms of H also reveals the property of G, that G = G~^. 
Finally, one may use the property that two square matrices related by a similarity 
transformation A' = SAS^^ obey det(/ + A') = det(i + A) to regroup the matrix 
products of (3.8) in a different manner if desired. 

The fact that the exponents in A^2 contain N+1 terms comphcates the derivation 
of a determinant formula for Z-i- 

3.3 Asymptotic behavior 

In this section we will examine how our solution behaves for x — > ±oo and t — > ±oo 
respectively. Since the dependence of our solutions on the worldsheet coordinates is 
encoded in the factors liZi^ the asymptotic behavior of the A^-magnon solution will be 
determined by their respective real parts. 
Using notation similar to [18], we define 

ui = i{Zi - Zi) = Kix - i/it, 
wi = Zi + Zi, 

vi = wi- t, (3.11) 



with 



XXi-Xi){l + \Xi\^) 2(l + rf)nsinf 



\1-Xif\l + Xif 1 + rf -2rf COS Pi' 
_ -^(A; - A^') _ 2nsinp; 
' |l-A,ni + A,p l + rf-2rfcosp,' ^""-'^^ 



^The matrix A is not to be confused with the Lagrange multipher of (2.1). 
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and in the second equality we have also employed the usual parametrization A; = ne^'^'-l'^ 
for the spectral parameters. Additionally, the relations (3.12) imply 

2iZi = ui + iwi, 2iZi = —ui + iwi- (3.13) 

The parameter range for a single dyonic magnon is r e (0, oo) and p e [0, 27r), with 
p ~ p + 27r for any other p. We can use the same restrictions for our parameters ri,pi 
of the A^-magnon solution, in which case the k,i are clearly positive. From the formulas 
(3.5)-(3.6) after wc factor out (3.7), it is then easy to see that the our solution has its 
boundaries on the equator of on the Zi plane. Namely, for x — > ±oo the boundary 
conditions (2.5) are satisfied, with p = Y^^=iPi ^ expected. 

Next, we proceed to determine the behavior of the solution for t — > ±oo and large 
magnon separation. Without loss of generality, we can assume that the magnons are 
ordered such that their velocities ^ obey 

^ > ^ > ... > ^. (3.14) 

In order to focus on the /c-th magnon, we keep Uk fixed as t — > ±oo. This means 
that X should scale bs x — ^t-\-'^ and in total the ui will behave as 

ui = ni[ \t + Hi — . (3.15) 

In particular, the limit i — > — oo under the aforementioned ordering and scaling implies 

Ui,U2, . . .,Uk-l +00, 

Uk finite, (3.16) 
Uk+i, Uk+2, . . . ,un ^ -oo. 

Thus, it is easy to see from (3.2)-(3.4) that the terms which dominate in the limit have 
fii = 1 for i G {1, . . . ,k — 1, k, k + 1, . . . , N} and i E {l,...,k — l,k,k + l,..., N} in 
the case of Ni, D, and i E {1, . . . , k — 1, k, k, k + 1, . . . , N} in the case of N2, with the 
rest of the /^.'s being zero. 

Up to common factors that will eventually cancel out (including the divergent 
terms) , we can express the limiting values of A^i , N2 and D as 

fc-l N 

Ni-^HhU'Xi {Xk f+ e"'= + ~Xk f- 6-"'=) e^"'^ (3.17) 



1=1 i=k+i 

k-l N 



N2-^Yl XiYl Xi A22/ie2™S 



1=1 i=k+i 
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where /+, /_, h are functions of the spectral parameters given by 

fc-l TV 



/+ = JJ^ I'^fc ~ '^/T Y\. i''^'^ ~ -^'1^' 

1=1 l=k+l 
fc-l N 



(3.18) 



1=1 
fc-l 



;=fe+i 



N 



h - n(^^ - ^i^Ch - A;) n (^fc - ^')(Afc - A,). 



1=1 



l=k+l 



Noticing that = f+f-, and with the help of (3.1), (3.17) and (3.18), we can write 
the i — > — oo limit of the A'^-magnon solution as 



Jt 



Z2 = e'^'e'"' 



cos ^ + i sin ^ tanh(iife + 5u^{k)) 



sm 



cosh [lifc + 5w_(A;)] ' 



where^ 



with 



^ ^ fc-l AT 

bu^{k) = - log ^ ^ ^Uk,l - ^ 5xife,/ 



buk^i = log 



i=fc+i 



Afc — A/ 



Afc — A^ 



(3.19) 



(3.20) 



(3.21) 



and the phase factors e*^^, e*^^ are independent of x and t. For completeness, we can 
write them explicitly as 



fc-l 



Vi[\\ n [t 



\\\ 



1/2 N 



;=fc+i 



A, 



A, 



1/2 



exp 



'fe-i 



AT 



9 E^^'- E 



./=1 



=fe+l 



(3.22) 



Equation (3.19) is precisely the single magnon solution on M x S*^ [8, 18], up to a pure 
phase and a shift in ti^, which reflects the additional freedom of the solution. 



"^The signs of 8u±{k) are chosen for compatibility with the most standard method of determining 
time delays, whereby one performs the ansatz = —i'kSt±(k) and solves for the position of the 
magnon's peak, given by —Vk6t±{k) + 5u±{k) — 0. Note the agreement with the definition (3.26) 
below. 
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The case t ^ oo can be treated in a similar manner, yielding (3.19) with 

5u_{k) — >• 5uj^{k) = —Su_{k), 
9i — > —9i, 



(3.23) 



(3.24) 



Since k is arbitrary, we have in fact proven that asymptotically our A^-magnon solution 
splits into N single magnon solutions. Each magnon retains its shape after scattering 
with the rest of the magnons, with the effect of the interaction being encoded only in 
a relative shift in Uk, 



5u{k) = 5u+{k) — 5u-{k) — —25u-{k). 
Because of (3.12), the shift in Uk is usually interpreted as a time delay [26], 

fc-l N 

- ^ (>tk,i + ^ Stk,i, 

1=1 l=k+l 

where 



St 



k,l 



= 2V- 



^k\ 1 1 + -^A 



log 



Afc — A/ 



Afc — A/ 



(3.25) 



(3.26) 



(3.27) 



is the time delay that occurs because of the interaction of the /c-th with the /-th magnon, 
namely two-magnon scattering. 

Hence, our A^-magnon solution exhibits the property of factorized scattering, as 
expected by the integrability of the R x 5"^ cr-model and its classical equivalence to the 
complex sine-Gordon system. Finally, the dyonic two-magnon time-delay we retrieved 
in (3.27) is in complete agreement with [27, 28, 19]. 
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A. Construction rules - Examples 

To help clarify the meaning of the formulas (3.1)-(3.4) and (3.5)-(3.6), we reduce them 
to a simple set of rules for the construction of A^i, A^2, D. These rules may also facihtate 
computer code for generating A^-magnon solutions. 
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The iV-magnon solution can be written as 



and it contains N spectral parameters Aj along with their conjugates Aj that we can 




arrange as the set A = {Ai, Ai, A2, . . . , A^v, \n}- 

In order to write the denominator D we take all the possible subsets of numbers 
of the set A. There are (2iV)!/iV!^ such subsets. For each subset we form a product 
and then D is the sum of all those products. Let us see how to form the product for a 
specific subset B. The product contains 

a) an exponential with exponent 2i - Z{Xi) = 2i Zi, where Aj are all the A's 
that belong to B, 

b) all the possible differences Aj — Aj, i < j, where Aj, Aj all belong to the subset 
B and 

c) finally all the possible differences Aj — Aj, i < j, where Aj, Aj all belong to the 
complement subset of B. 

The rules for Ni are the same as D except that now the product contains in addition 
all the A's that belong to the subset B. 

The rules for N2 are the same as the rules for A^i, but now all the subsets B should 
have + 1 elements instead of N and the product contains all the A's that belong to 
the complement subset of B instead of the B itself. 

As an example let us write A^i, N2, D in the case of 1, 2 and 3-magnons. 
For 1-magnon we have [8] 



= Aie^'^i + Aie^^^i 
N2 = Xue^'^^'+^'\ 



(A.2) 



For 2-magnons we have [18] 



D = AnA22e'^(^^+^^) + Ai2Ai2e^^(^^+^^) + Ai2Ai2e2*(^^+^^) 
+Ai2Ai2e2^(^^+^^) + Ai2Ai2e^^(^i+^^) + A22Ane^^(^^+^^), 



A^i = AiAiAnA22e'^(^^+^^) + AiA2Ai2Ai2e2^(^i+^^) + AiA2Ai2Ai2e2^(^i+^^) 

+ AiA2Ai2Ai2e''(^^+^^) + AiA2Ai2Ai2e''(^i+^^) + A2A2A22Ane'^(^^+^^), 



N2 = A2AnAi2Ai2e2^(^i+^i+^^) + A2AnAi2Ai2e^^(^i+^i+^^) 
+AiAi2Ai2A22e'^(^^+^^+^^) + AiAi2Ai2A22e^^(^^+^^+^^) 



(A.3) 
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For 3-magnons we have 

D — AiiAi2Ai2A23A23A33e^*''^^'''^^^^^^ + AiiAi2Ai2A23A23A33e^**-^^'''^^'''^^'' 

+ AiiAi3Ai3A22A23A23e^**-^^'''^^^^^^ + AiiAi3Ai3A22A23A23e^*''^^'''^^^^^^ 
+ Ai2Ai2A22Ai3Ai3A33e^'*-'^^'''^^'''^^^ + Ai2Ai3A23Aj2Aj3A23e^*^^^'''^^'''^^^ 
+ Ai2Ai3A23Ai2Ai3A23e^*^'^^''"^^''"^*^ + Ai2Ai3A23Ai2Ai3A23e^*'^'^^''"'^^''"'^^^ 

+ Ai2Ai3A23Ai2Ai3A23e^*^'^^^^^^^^^ + Ai3Ai3A33Ai2Ai2A22e^*^'^^^'^^^'^^^ 
+ Ai2Ai2A22Ai3Ai3A33e^**-^^'^'^^''"'^^^ + Ai2Ai3A23Ai2Ai3A23e^*''^^'^'^^"^'^^^ 
+ Al2Ai3A23Ai2Ai3A23e^**-^^"'"'^^"'"^^^ + Ai2Ai3A23Ai2Ai3A23e^*''^^'*"'^^'*"'^^^ 
+ Ai2Ai3A23Ai2Ai3A23e^**'^^'''^^'''^^'' + Ai3Ai3A33Ai2Ai2A22e^'*-^^''''^^''''^^'' 
+ A22A23A23AllAi3Aj3e^**-^^^^^'^^^^ + A22A23A23AiiAi3Ai3e^*^^^'''^^'''^^^ 
+ A23A23A33AiiAi2Ai2e^'*-^^'''^^'''^^^ + A23A23A33AiiAi2Ai2e^*^^^'''^^'''^^\ 

Ni — AiAiA2AijAi2Aj2A23A23A33e^*'-^^'''^^'''^^^ + AiAiA2AixAi2Ai2A23A23A33e^**-^^'''^^'''^^^ 
+ AiAiA3AiiAi3Ai3A22A23A23e^**-^^'''^^^^^'' + AiAiA3AiiAi3Ai3A22A23A23e^**'^^^^^^^'^'' 

+ AiA2A2Ai2Ai2A22Ai3Ai3A33e^'^'^^'''^^'''^^^ + AiA2A3Ai2Ai3A23Ai2Ai3A23e^'^'^^'''^^'''^^^ 
+ AiA2A3Ai2Ai3A23Ai2Ai3A23e^*^'^^'^'^^^'^^^ + AiA2A3Ai2Ai3A23Ai2Ai3A23e^*^'^^"^'^^'^'^^^ 
+ AiA2A3Ai2Ai3A23Ai2Ai3A23e^*^'^^''"'^^''"'^^^ + AiA3A3Ai3Ai3A33Ai2Ai2A22e^*^'^^"'^'^^^'^^^ 
+ AiA2A2Ai2Ai2A22Ai3Ai3A33e^*^'^^'''^^"'"^^^ + AiA2A3Ai2Ai3A23Ai2Ai3A23e^*^'^^"'"'^^''"'^^^ 

+ AiA2A3Ai2Ai3A23Ai2Ai3A23e^**-'^^'*''^^'''^^'' + AiA2A3Ai2Ai3A23Ai2Ai3A23e^'*-'^^^'^^"'"'^^^ 

+ AiA2A3Ai2Ai3A23Ai2Ai3A23e^**''^^^^^'''^^'' + AiA3A3Ai3Ai3A33Ai2Ai2A22e^**'^^^'^^"'"'^^^ 
+ A2A2A3A22A23A23AllAi3Ai3e^**-^^'''^^'''^^'' + A2A2A3A22A23A23AiiAi3Ai3e^**-^^'''^^'''^^^ 
+ A2A3A3A23A23A33AiiAi2Ai2e^'*^^^'''^^'''^^^ + A2A3A3A23A23A33AiiAi2Ai2e^'*^^^'''^^'''^^\ 

N2 = A3A3AiiAi2Ai2Ai2Al2A22A33e2*(^i+^i+^=+^^) + A2A3AiiAi2Ai3Ai2Al3A23A23e2*(^i+^i+^=+^3) 

+ A2A3AiiAi2Ai3Ai2Ai3A23A23e2^(^i+^i+^^+^^) + A2A3AnAi2Ai3Ai2Ai3A23A23e2^(^^ 

+ A2A3AiiAi2Ai3Ai2Ai3A23A23e^*^'^^~'~'^^~'~'^^~'''^^^ + A2A2A11 Ai3Ai3Ai3Ai3A33A22e^*^'^^~'~'^^~'~'^^~'~'^^' 
+ AiA3Ai2Ai2Ai3A22A23A23Ai3e2'(^i+^^+^^+^3) + ^^^g^^^^^_^^_^^_^^_^_g^_^g2i(Zi+2.+2.+23) 

+ AiA2Ai2Ai3Ai3A23A23A33Ai2e2'(^i+^^+^3+23) + AiA2Ai2Ai3Ai3A23A23A33Al2e2*(^i+^=+^='+^=') 

+ AiA3Al2Ai2Ai3A22A23A23Ai3e2'(^i+^^+^^+^3) + AiA3Ai2Ai2Al3A22A23A23Ai3e2*(^i+^=+^^+^=') 

+ AiA2Ai2Ai3Ai3A23A23A33Ai2e''(^^+^^+^^+^^) + AiA2Ai2Al3Ai3A23A23A33Ai2e^*(^i+^=+^3+^=') 

+ AiAiA22A23A23A23A23A33Aiie^*''^^~'~'^^~'~'^''~'''^^-', 

(A.4) 
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where Xij = Aj — Xj, and Zi = zj (Aj — 1) + zjiXi + 1). 
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